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We have identified anomalous behavior of the escape rate out of the zero-voltage state in Joseph¬ 
son junctions with a high critical current density J c . For this study we have employed YBa 2 Cu 3 07 _a; 
grain boundary junctions, which span a wide range of J c and have appropriate electro-dynamical 
parameters. Such high J c junctions, when hysteretic, do not switch from the superconducting to 
the normal state following the expected stochastic Josephson distribution, despite having standard 
Josephson properties such as a Fraunhofer magnetic field pattern. The switching current distribu¬ 
tions (SCDs) are consistent with non-equilibrium dynamics taking place on a local rather than a 
global scale. This means that macroscopic quantum phenomena seem to be practically unattainable 
for high J c junctions. We argue that SCDs are an accurate means to measure non-equilibrium 
effects. This transition from global to local dynamics is of relevance for all kinds of weak links, 
including the emergent family of nano-hybrid Josephson junctions. Therefore caution should be 
applied in the use of such junctions in, for instance, the search for Majorana fermions. 

PACS numbers: 74.50.+r, 85.25.Cp, 74.40.Gh, 74.78.Na 


I. INTRODUCTION 

Rich physics, such as the Josephson effect 1 } quantum 
coherence® and quantum interference, arise when two 
coherent quantum systems are weakly coupled. Some ex¬ 
amples of systems that exhibit such physics are supercon¬ 
ducting Josephson junction^, the flow of superfluid 4 He 
through an array of nano-aperture^, and the observation 
of the Josephson effect in Bose-Einstein condensate^. 
Superconductors have traditionally held a special place 
as test-bench systems for these kinds of quantum phe¬ 
nomena due to the ease of scaling and integrating them 
into real quantum devices, offering a high degree of mea¬ 
surability and control of a macroscopic wave function. 

Every experiment or application using a superconduct¬ 
ing weak link is based on how the phase differe nce <p be¬ 
tween the electrodes evolves in time and spacd^yyiU The 
large variety of barriers now available between the su¬ 
perconducting electrodes offer novel functionalities and 
efficient tuning of physical processes occurring at the 
nanoscale and at different interface^®!] A recent ex¬ 
ample is the proposal to use the Josephson effect for the 
detection of the Majorana fermions 14 . 

Progress in material science in producing a larger 
variety of interfaces and in nanotechnologies applied 
to superconductivity, is promoting a rethinking of the 
phase dynamics of Josephson junctions (JJs). Here we 
give evidence of a breakdown of a fundamental tenet 
of the Josephson effect: the transition from the super¬ 
conducting to the normal state does not follow the ex¬ 


pected stochastic Josephson phase dynamics, but has a 
more intriguing balance between local and global en¬ 
ergy processes. The onset of non-equilibrium effects 
is the key to describe the local processes, which may 
radically change perspectives on how to interpret the 
physics and make predictions on the performances of a 
large variety of ” smart” Josephson devices. We mea¬ 
sure switching current distributions (SCDs) which cod¬ 
ify the very general process of the escape of a particle 
(phase) from a potential well in a JJ 2 ^, keeping track 
in our case of non-equilibrium effects. Roughly speak¬ 
ing, SCDs are obtained in JJs with hysteretic current- 
voltage (I-V) characteristics by counting the number of 
times the system switches from the superconducting to 
the resistive state within a small window of bias cur¬ 
rent, when ramping forth and back the bias current. 
Thermally activated processes are well understood in JJs 
both in the underdamped 15 17 and in the moderately 
dampecP^ regime. The transition to the macroscopic 
quantum tunneling regime has been theoretically®* 21 
and experimentally 22 2 ^ widely investigated. SCD mea¬ 
surements focus on ”the very moment” at which resis¬ 
tance originates in superconducting weak links. Thus we 
use the power of encoding information of non-equilibrium 
local processes in fluctuations to characterize these dy¬ 
namical processes. 

In the resistively and capacitively shunted junction 
(RCSJ) modeP, the damping parameter Q = uj p RC is 
proportional to the square root of I c via the plasma fre¬ 
quency uj p = (2eI c /hC) 1 ^ 2 at zero bias current, where 
R and C are the resistance and capacitance, respec- 
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FIG. 1: a) Three-dimensional representation of the SCDs measured for various temperatures on three GB YBCO JJs. b) At 
low J c values (5-10 2 A/cm 2 ), the SCDs of JCT A are confined to a small range of currents, with a of the order of 10 nA. With 
increasing J c the histograms progressively cover a larger interval of currents, both absolutely and relatively (JCTs B and C in 
panels c and d respectively). In panel e) the escape rate curves r(7) computed from the SCDs of JCT C have been shown for 
the same temperature range as in panel d. 


tively. In a mor e general approach, Q has a frequency 
dependenc e] 18 1 25 1 , which includes the effects of the exter¬ 
nal shunting impedance. A junction cannot sustain an 
unlimited increase in the critical current I c and thus 
in the quality factor Q through larger critical current 
density J c while still preserving all the properties of the 
Josephson effect and all the features of the underdamped 
regime in the I-V curves. 

The breakdown of the Josephson dynamics is found at 
higher values of the critical current density. The clas¬ 
sical Josephson phase dynamics, which takes place in 
junctions characterized by lower critical current densi¬ 
ties J c , are replaced at high J c values by a regime driven 
by non-equilibrium dynamics where phase information 
is lost. Non-equilibrium effects produce hysteretic I-V 
characteristics and modify the influence of dissipation, 
thus becoming measurable through modeling of the SCD 
in terms of heating modes. The heating modes follow 
a heat diffusion-type equation in analogy to phase slip 
eventJ- 26 27 . In the present work, this transition from 
classical to non-equilibrium phase dynamics is found for 


high critical temperature superconductors (HTS) grain 
boundary (GB) junctions, but should be expected for any 
kind of Jjd^C 29 . Specific thresholds may depend on the 
type of junctions and materials, but the features of the 
transition are universal. HTS GB junctions are the ideal 
system to identify this transition because of the possibil¬ 
ity of varying the c ritical supercurrent density J c over a 
wide range of valueP^M 


II. JOSEPHSON PHASE DYNAMICS 

YBa 2 Cu 307 _ iC (YBCO) off axis biepitaxial GB junc¬ 
tions provide a large variety of transport regimes because 
of d-wave order parameter effects as well as the modula¬ 
tion of the barrier transparency through different relative 
orientation of the electrodes 30 32 . Details on the prepa¬ 
ration and properties of t he GB junctions used in this 
work are given in RefsP-^^IS] p resen t here data from 
junctions (JCT) with low (JCT A), intermediate (JCT 
B) and high (JCT C) values of J c . For each junction J c 
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FIG. 2: a) Temperature dependence of a measured on GB JCTs A (blue squares), B (red triangles) and C (black points) 
respectively. The blue and red solid lines are Monte Carlo simulations of the phase dynamics, according to multiple escape and 
retrapping processes in the washboard potential (shown in panel b), with Q = 1.30 and 1.17 for JCTs A and B respectively. 
Some data from sample A have been presented in Ref. 40 previously, c) Simulations of the thermal dependence of cr as a 
function of J c for different values of the Q damping parameter (full color lines), confined to the moderately damped regime, are 
compared with experimental data of JCT C (black points). Keeping all the other junction parameters fixed, an enhancement 
of J c leads to an increase of I c and Q. An increase in Q produces steeper a tails above T* and cannot reproduce the broadened 
experimental data of JCT C. This is even more evident in the two dimensional (cr — T) projection. For completeness, in the 
(cr — T) projection the color points referred to the results of Monte Carlo simulations have been reported for Q values ranging 
from 1 to 10, and compared with the data of JCT C (black points). The colored lines are guides for the eye. At high J c the 
Langevin approach does not hold anymore and non-equilibrium concepts apply. Simulations fitting the SCDs of JCT C (see 
Fig.@ are consistent with multiple heating events as reported pictorially in panel d. As explained in the text, 77 is the bath 
temperature and Tth is the threshold temperature above which the transition to the resistive state occurs. In panel d) the 
arrows qualitatively sketch the temperature jump due to a single heating event and the number of events necessary to induce 
the switching to the resistive state. 


has been determined as the ratio between the measured 
I c in the I-V curves and the geometrical area of the junc¬ 
tion (width x thickness of the YBCO film)^. For JCT C 
we have employed nan ofabri cation techniques to be in the 
appropriate range of /j 25 | 35 i To study their escape rates 
we have thermally anchored the samples to the mixing 
chamber of a He 3 /He 4 Oxford dilution refrigerator and 
performed measurements of SCDs. Measurement proce¬ 
dures and filtering used in this experiment are described 
elsewherJSO. 

Switching histograms of JCTs A, B and C are shown 
in Fig. [l] The 3-D view gives an intuitive picture: SCDs 


cover distinct current ranges, and when J c increases, they 
become broader. In the right part of the figure, each 
set of SCDs is displayed with appropriate scales for a 
better view of the details of their temperature depen¬ 
dence. In Fig. the escape rate V out of the zero- 
voltage state as a function of the bias current / has been 
plotted for JCT C in the same temperature range. T(I) 
curves have been computed from the SCDs following Ful¬ 
ton and Dunklebergei^ 5 . The standard deviation a of the 
histograms is reported for all junctions in Fig. |2^i. The 
combined analysis of cr, of the skewness 7 = m 3 /cr 3 (m 3 
being the third central moment of the distribution), and 
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of their temperature evolution characterize the phase dy¬ 
namics. 

The histograms from JCTs A and B match well the 
predictions of the RCSJ model for moderately damped 
JJs. Low values of especially in HTS d- 

wave junctions, where intrinsic sources contribute to 
dissipatio n 2 - 5 -!— I, lead to the moderately damped (MD) 
regime (Q ~ 1). In the MD regime the SCDs become 
narrower with increasing temperature for temperatures 
above T*, which is defined as the temperature at which 
cr(T) is largest (in the case of JCT B T* is about 1.0 K). 
This is often repo rted in the literature as the phase dif¬ 
fusion (PD) regimd^ ^ 38 l 4Q l , indicating a phase diffusive 
dynamics due to multiple escape and retrapping events in 
the washboard potential (sketched in Fig. or equiv¬ 
alently described in terms of thermal activation above 
a dissipation barrier 39 . Below T*, SCDs obey the ther¬ 
mal activation (TA) regime. Upon further lowering the 
temperature the junction can event ually enter the macro¬ 
scopic quantum tunneling regime 19 21 (see Fig. [TJd, 
JCT A). The MD regime gives a very important refer¬ 
ence for our discussion of sample C. In both cases the 
width of the SCD decreases upon increasing the temper¬ 
ature, as shown in Figs. [l]and[2] Only an exact fitting of 
the SCDs at different temperatures can reveal the very 
different physics occurring in the various junctions. 

JCTs A and B of the present work have been chosen 
since they are fabricated on different substrates (see Ref. 
[33) and both show the transition to the PD regime, in 
which the derivative da/dT is negative, but with a differ¬ 
ent detailed T dependence 36 ! The choice of the substrate 
affects the effective capacitance of the circuit the junc¬ 
tion is embedded into 41 . As a consequence of this choice, 
for JCT B we observed the transition from TA to the PD 
regime, while for JCT A the TA regime is completely sup¬ 
pressed and we observed a direct transition from MQT 
to the PD regime 40 . We show below that the anomalous 
thermal behavior of high J c JCT C cannot be explained 
in terms of frequency dependent damping or as a conse¬ 
quence of the shell circuit, and the numerical simulations 
of Fig. fully support this. 

In JCTs A and B the rate of decrease in a above T* 
with increasing T is well described by the Monte Carlo fit 
of the phase dynamics® (blue and red solid lines in Fig. 

respectively). The temperature dependence of the 
skewness 7 gives additional distinc tive cr iteria to confirm 
the intermediate dissipation regim^®®. This is signaled 
by a transition for JCT B for instance from 7 -0.9 at 

low temperatures (T ~ 0.1 K) in the TA to 7 ^ -0.1 above 
the transition temperature T* ~ 1.0 K in the PD, indi¬ 
cating a progressive symmetrization of the SCDs® (see 
Fig. §. Therefore, the moderately damped regime gives 
the opportunity of comparing cr(T) and y(T) dependen¬ 
cies of SCD spectra at higher temperatures, which turn 
out as unambiguous distinctive criteria for the different 
switching modes. 

JCT C is characterized b y hig h values of J c close to 
those observed in nanowires 43 44 . This device exhibits 
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FIG. 3: Temperature dependence of the skewness 7 of the 
SCDs measured on JCTs B (red triangles) and C (black cir¬ 
cles) respectively. Retrapping processes in the PD regime 
cause a progressive symmetrization of the SCDs of JCT B, 
signaled by the thermal dependence of 7 , which is almost 0 
above the transition temperature T*. Asymmetric SCDs, as 
the histogram measured at T = 0.1 A < T* (shown in the top 
left corner of the figure), become Gaussian-like distributions 
above T* (see the histogram measured at T — 4 K > T* in 
the top right corner). The histograms of JCT C are asymmet¬ 
ric over the whole temperature range, including for temper¬ 
atures above T** (as shown by the SCD measured at T = 7 
K reported in the lower right corner of the figure). The 
black dashed line indicates the mean value of the skewness 
(7--0.8). 


radically different phase dynamics above 3 K, which rep¬ 
resents a transition temperature we indicate as T**. Also 
in the case of JCT C, T** is defined as the maximum of 
cr(T), but we prefer to use a different symbol to stress the 
distinct switching dynamics, as pointed out below. The 
rate of decrease of a above T** turns to be a distinctive 
marker of the phase dynamics. In device C the slope of 
cr(T) above T** is much smaller when compared to those 
of JCTs A^ni and B and of all moderately damped JJs. 
The smooth decrease of a for T > T** cannot be de¬ 
scribed in terms of the intermediate dissipation regime, 
as evident from the numerical simulations reported in 
Fig. According to the RCSJ model, keeping all the 
other junctions parameters fixed, an increase of J c leads 
to an enhancement of I c and Q. The increase of Q moves 
T* to higher values, and the negative slope of cr(T) above 
T* becomes steeper and steeper (as shown by the numer¬ 
ical simulations for Q ranging between 1 and 10 reported 
in Fig. [2J3). 

In addition, while moderately damped JJs show a pro¬ 
gressive symmetrization of the switching histograms near 
to and above T* (see JCT B in Fig. |3]), SCDs of JCT 
C are asymmetric over the entire temperature range. In 
JCT C, the 7 factor is temperature independent, consis¬ 
tent with what is observed in pure phase slip systems 42 
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(see Fig. [3| the black dashed line indicates the mean 
value 7 c±-0.8 for JCT C). These behaviors are quite dis¬ 
tinct and do not fall in the framework of any regime of 
the RCSJ model. 


III. SWITCHING DYNAMICS OF HIGH 

CRITICAL CURRENT DENSITY JOSEPHSON 
JUNCTIONS 

We find that the numerical simulation of a transition 
driven by local heating events accounts well for devices 
in the J c interval ( 10 4 A/cm 2 - 10 5 A/cm 2 ), as s hown for 
JCT C in Fig. if- Hysteresis in I-V curvej ^ 45 * 46 * does 
not necessarily indicate canonical Josephson phase dy¬ 
namics, even in the presence of a Fraunhofer magnetic 
field pattern 29 . It may rather arise as a result of local 
heating processes, possibly induced by intrinsic inhomo¬ 
geneous composition unavoidable for high J c junctions. 
Large values of Joule power density deposited in the weak 
link can induce a self-heating process during the switch 
to the resistive branclF^. The absence of a set of self- 
consistent electrodynamics parameters to describe JCT 
C is a strong indication of the failure of the standard 
Josephson dynamics. This failure is of ge neral relevance, 
applying both to conventional low T c JJ d — * and to the 
emergent class of hybrid nanoscale junctions 10 12 JQ2[ For 
larger values of J c , heating driven mechanisms become 
dominant with a transition to the normal state locally 
in the junction area. These events are the mirror of 
non-equilibrium processes and can be modeled as ” phase 
slips entities” (PSEs), as confirmed by the details of the 
simulations which are reported in Refs. 133I47II491 in 
the sense that they are local processes, break the co¬ 
herence of the phase information and are described by 
a heat diffusion-like equation. In particular, the prob¬ 
ability for a single heating event can be still described 
in terms of the Langer-Ambegaokar-McCumber-Halperin 
(LAMH) theor}!^ and further extension^!, and PSEs 
can still be approximately assumed to be far apart in 
time. 

The typical time scale of a PSE is of the order of the 
Ginzburg Landau relaxation time ~ r(u) = (1 — u)~ x tq 
where To = irh/8kBT c and u = T/T c < 1, with T c the 
superconducting critical temperature. The energy dissi¬ 
pated by a PSE Epse{I ) = <\> 0 d is rather high, as the 
modulus of the order parameter has to vanish over a 
length £(T) = £ 0 (1 — Ti 4 ) 1 / 2 /(1 — u 2 ), where £o is the zero 
temperature coherence length in the superconductor and 
(fio is the flux quantum. Following an approach proposed 
in Ref. [47] for low critical temperature superconducting 
(LTS) wires, our numerical simulation of the tempera¬ 
ture jump induced by a PSE obeys the phenomenological 
diffusive equation for the relaxation of the temperature 
gradient: 

d E+a{T,T b ) ST = r(T b ,t) + V (T,I)J2Ht-U). (1) 


Here ST = T — T& is the deviation from the tempera¬ 
ture of the bath T^. The relaxation coefficient (T, T&) 
depends on the thermal conductivity K{T ), on the ther¬ 
mal capacity C V (T ) of the phase slip volume, and on X 5 33 . 
r(T^,t) is the noise source due to the environment with 
an admittance Y(cj), while 7 ?(T, I) is the temperature 
jump due to the PSEs which occur at the stocasticaly 
distributed times tj® 3 . ^(T, I) is implicitly defined by 

the following equation: 



In Fig. [4^i we show the SCDs derived from the ex¬ 
periment on JCT C over a wide temperature range, 
T E [0.25A, 10K]. The continuous red curves corre¬ 
spond to the fit obtained integrating eq. [l] To quantify 
how many heating events are needed to escape to the fi¬ 
nite voltage state, we define the threshold temperature 
T t h , i.e. the temperature above which the system is defi¬ 
nitely out of the zero voltage state. At low temperatures, 
T < T**, our fitting procedure does not depend on T t ^ . 
This is the regime when a single heating event is enough 
to drive the transition. Following Ref. 071 we assume 
that both the specific heat and the thermal conductiv¬ 
ity K(T ) are weighted averages of their BCS and Fermi 
liquid limits. Two main effects discriminate between the 
low temperature and the high temperature behavior. 

The temperature jump 77 (F) depends on temperature 
because the specific heat is strongly temperature depen¬ 
dent. At low temperatures the specific heat is quite low, 
thus with each PSE there is a considerable increase in the 
temperature. In addition, the thermal conductivity and 
the thermal capacity are both quite low as the system is 
deeply into the superconducting phase. In a rough sense, 
the system is rather isolated from the environment and 
the thermal shock due to a heating event is destructive 
for the superconducting state. For this reason a single 
heating event can induce a direct jump to the resistive 
state at low temperatures: it induces a large local heat¬ 
ing that is difficult to dissipate. The system is not at 
equilibrium with its environment, and we can define an 
effective temperature T eff_ for the system, which is higher 
than Tb as shown in Fig. &• 

At high temperatures we are in the opposite regime 
of small 77 (T) per heating event 33 . In addition both 
the thermal conductivity K(T ) and the thermal capac¬ 
ity C V (T ) increase with increasing temperature as well. 
Thermal diffusion and contact with the environment is 
more effective and multiple PSEs are required for switch¬ 
ing. This occurs above T**, where the derivative dcr/dT 
is negative. T e ff and T^ coincide above T** (see Fig. 
[4|:), which we interpret as the temperature at which the 
system is able to thermalize during the time interval be¬ 
tween well separated heating events. In this temperature 
range, the value chosen for the fitting parameter T t h be¬ 
comes important to quantify the number of successive 
PSEs, which are responsible for the transition. In Fig. 
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FIG. 4: a) SCDs measured on JCT C (black circles) along 
with the fit (red solid lines) resulting from numerical simu¬ 
lations of heat diffusion-like dynamics, according to eq. [I] 
b) Fitting procedure above T**: single heating event dynam¬ 
ics does not reproduce the experimental histogram, signaling 
a multiple PSE process. Dissipation due to heat relaxation 
plays a relevant role, c) Temperature dependence of T e // es¬ 
timated from numerical simulations of eq. 1. Below T**, a 
single heating event induces the switching and T e // > TJ,. 
Above T**, T e ff and Tb coincide since the system is able to 
thermalize in the time interval between well separated heating 
events, and multiple PSEs are needed to overcome T t h- 


[4]3 the difference between a single and multiple heating 
events is shown. At 7 K for instance about 4-6 PSEs 
are needed to reproduce the experimental SCD. Exper¬ 
imental SCDs complemented by the numerical simula¬ 
tions follow the passage from single to multiple heating 
events. T his is si milar to what has been observed in LTS 
nanowire d 43 * 44 * 47 . A consistent set of the junction param¬ 
eters (temperature jump 77, number of heating events) 
can be extracted from these simulations, as discussed in 
Ref. [33l 


IV. DISCUSSION AND CONCLUDING 
REMARKS 

The passage above a characteristic J c threshold to a 
state that does not sustain Josephson phase-coherence 
should be universally expected in weak link systems. The 
threshold probably depends on the type of junction, re¬ 
sponding to structurally different ways of enhancing bar¬ 
rier transparency and generating sources of local heating. 
The threshold of 10 4 — 10 5 A/cm 2 found for our j unc¬ 
tions jsconsistent with the general empiric rule on HTS 
GB d^ —! In bicrystal junctions, an increase in the mis- 
orientation angle between the two electrodes determines 
a well-known decrease of the critical current density. For 
small misorientation angles devices are found in the flux- 
flow regime. For 6 > 20° Josephson phenomenology 
starts to appear for J c values of 10 3 A/cm 2 ^ 31 ! Several 
effects have also been observed in conventional trilayer 
Nb/Al-Oxide/Nb JJs above the same threshold of about 
10 4 A/cm 2 that could be attributed to non-equilibrium 
dynamic d 28 l 29 [ 

Some analogies can be established with what is ob¬ 
served in 4 He superfluid 6 1 Here the passage from weak to 
strong coupling manifests itself in a change in the current- 
phase relation [I(p)]. In the strong coupling regime, 
where the healing length of the superfluid is lower than 
the diameter of the constriction, deviations from the sin p 
relation appear, while sinusoidal Josephson oscillations 
have been measured in the opposite limit^. An increase 
in J c and in the coupling between the electrodes leads to 
the presence of other harmonics in the I(p), which might 
become multivalued 2 ! 

A heat diffusion-like model breaking phase-coherent in¬ 
formation is consistent with our data on high J c JCT C. 
Here the switch to the normal state is accompanied by a 
local release of energy characteristic of a PSE. When de¬ 
parting from the supercurrent branch, non-equilibrium 
processes produce an unexpected heating. This is sur¬ 
prisingly different from what commonly is accepted for 
hysteretic Josephson junctions, where heating can only 
influence the retrapping phenomena, a s a m emory of the 
history of heating in the resistive statd® 51 . 

In the framework of the RCSJ model the resistance 
arises from non-local properties well described by the 
(frequency dependent) quality factor Q of the circuit. 
This reminds us of the Landauer picture of quantum 
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transport, in which quantum interference acts at the in¬ 
terface while the dissipation and memory loss occurs in 
the contacts. On the other hand, when the switch is 
driven by PSEs, it is the local properties of the order pa¬ 
rameter which are important: phase memory is lost when 
the modulus of the order parameter becomes zero. The 
model based on local heating events can be in principle 
further extended to extract the I (ip) and to fully define 
the role of dissipation in high J c junctions. 

To conclude, standard phase dynamics of a hysteretic 
Josephson junction collapses and cannot be sustained 
above some threshold of the critical current density J c . 
In high J c devices information is lost in non-equilibrium 
dynamics, and can be partly codified in the local heating 
process. This is of great relevance for all the experi¬ 
ments using low-dimensional barriers, which should be 
concerned about possible heating effects, leading to dis¬ 
torted phase information. Non-equilibrium effects would 
obviously invalidate a large number of key predictions for 
Majorana fermions and nanoscale superconductivity in 


Josephson junctions, since for instance they can give rise 
to zero bias anomalies in conductance measurements 52 . 
The proof of quantifiable non-equilibrium processes al¬ 
ready in the thermal regime poses severe constraints on 
the possible occurrence of macroscopic quantum phenom¬ 
ena at lower temperatures in high J c samples through 
standard SCD measurements. 
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